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The Sznajd model is a sociophysics model, that is used to model opinion propagation and consensus formation 
in societies. Its main feature is that its rules favour bigger groups of agreeing people. In a previous work, we 
generalized the bounded confidence rule in order to model biases and prejudices in discrete opinion models. In 
that work, we applied this modification to the Sznajd model and presented some preliminary results. The present 
work extends what we did in that paper. We present results linking many of the properties of the mean-field fixed 
points, with only a few qualitative aspects of the confidence rule (the biases and prejudices modelled), finding 
an interesting connection with graph theory problems. More precisely, we link the existence of fixed points with 
the notion of strongly connected graphs and the stability of fixed points with the problem of finding the maximal 
independent sets of a graph. We present some graph theory concepts, together with examples, and comparisons 
between the mean-field and simulations in Barabasi-Albert networks, followed by the main mathematical ideas 
and appendices with the rigorous proofs of our claims. We also show that there is no qualitative difference in 
the mean-field results if we require that a group of size q > 2, instead of a pair, of agreeing agents be formed 
before they attempt to convince other sites (for the mean-field, this would coincide with the q-voter model). 
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I. INTRODUCTION 

In the last years, the interest in interdisciplinary problems 
has increased among physicists, creating many research areas. 
One of these areas is sociophysics, that studies how assump- 
tions about the behaviour and social interactions of people in 
a "microcospic level" creates emerging social behaviours, like 
opinion propagation, consensus formation, properties of elec- 
tions, how wealth is distributed in society, among other top- 
ics. Typical approaches include modelling using deterministic 
celular automata, Monte Carlo simulations of models derived 
from ferromagnetic models (usualy Ising and Potts), mean- 
field approaches and diffusion-reaction processes 
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The Sznajd model is an opinion propagation model, orig- 
inally inspired by the Ising model in a linear chain, and is 
typically used to model consensus formation. It has spawned 
many variations, including the addition of noise, contrarian- 
like agents and undecided voters; as well as generalizations to 
more than 2 states (opinions) and to arbitrary networks [H U^ll . 
In all these variations, the most defining aspect of the Sznajd 
model is that it gives a greater convincing power to bigger 
groups of agreeing agents. Even though the importance of 
this effect has been known by psychologists since the 1950s 
Oh, it is often overlooked in other opinion propagation mod- 
els, for the sake of simplicity (this happens for example in the 
voter and in the Deffuant models ifl fllh ). 

In a recent work, we took the bounded confidence rule (that 
roughly says that people are only allowed to change opinions 
in a smooth way) that is common to many opinion propaga- 
tion models l^ flol[l5Tl . including the Sznajd model, and we 
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generalized it to model biases and prejudices in discrete opin- 
ion models (these generalized rules will be called by the um- 
brella term confidence rules). We applied this generalization 
to the Sznajd model and studied mainly the case with 3 opi- 
nions. In that work, we found a good qualitative (and in some 
cases quantitative) agreement between the model simulated in 
Barabasi-Albert (BA) networks fH and the mean-field equa- 
tions, but some of the results about the mean-field were still 
rather sketchy. In the present paper we give rigorous proofs 
(the main mathematical ideas are in a sepparate section and 
the detailed proofs can be found in the appendices) about the 
structure, existence and stability of the fixed points for the 
mean-field version of the Sznajd model with general confi- 
dence rules, finding a connection between these properties and 
graph theory problems using a graph derived only from qual- 
itative properties of the confidence rule. 

The results have some counterintuitive aspects and as such 
we provide both numerical solutions for the mean-field equa- 
tions and Monte Carlo simulations for the model in a BA net- 
work. In our calculations for the mean-field, we use a variant 
of the model, where at each timestep we choose q agents and if 
they agree, they attempt to convince r other agents; and show 
that there is no qualitative difference between all the cases 
with q > 2 (which includes the usual definition of the Sznajd 
model in an arbitrary network, with an arbitrary number of 
opinions, as we have considered in lfl7h ). 

n. THE SZNAJD MODEL WITH CONFIDENCE RULES 

The Sznajd model is an agent based sociophysics model for 
opinion propagation. In this model, a society is represented 
by a network (that is, a collection of nodes linked together by 
edges), where each node represents an agent (person), each 
edge is a social connection (friendship, marriage, acquain- 
tances, etc.) and each node i possesses an integer cr,-, between 
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1 and M, representing its opinion. In our generalization of the 
Sznajd model, as defined in lfl7ll . we introduce a set of param- 
eters Pa-^a-' (that are fixed and completely independent with 
the state of the network), and at each time step the following 
update rule is used: 

• A node i is chosen at random, and then a neighbour j 
of i is chosen. 

• If they disagree (cr^ ^ o~j), nothing happens. 

• If they agree, a neighbour k of j is chosen and is con- 
vinced of opinion o\ with probability p (Tfc _ >(Ti . 

We can interpret the first step as a conversation between two 
people that know each other, where they discuss some issue. 
If they disagree, none manages to convince the other. But, 
if they agree, they may set to convince another person that 
one of them knows and this person is convinced with a certain 
probability that depends only of its current point of view and 
of the opinion the pair is trying to impose. 

In the original model the probability weights p a ^a' are 
not dependent on a and a' . The reason why this probabil- 
ity should depend on both opinions is that, usualy an opin- 
ion includes prejudices about differing points of view (this is 
strongly related with the idea of cognitive dissonance in psy- 
chology [|7l |8l ll3IO . This generalization allows for complex 
interactions among the opinions in an unified way and can 
be seen as a generalization of the bounded confidence rules 
fiEa, as those rules can be recovered as special cases. Other 
modifications of the model can also be obtained this way: 

• When Pa^a' = 1 V a, a 1 we have the usual model. 

• If \a— u'\ < e => Pa^a' = 1 andp^o-' = otherwise, 
we have bounded confidence with threshold e. 

• Undecided agents can be modelled by a special state a, 
such that Pa'^a = V a' (undecided agents can only 
be convinced). 

• Cyclic interactions, like rock, paper, scissors (A con- 
vinces only B, that convinces only C, that convinces 
only A). 

This generalized version of the model has M{M — 1) pa- 
rameters, where AI is the number of opinions {p a ^a is irrel- 
evant and can always be taken as 0). These parameters can be 
thought as the elements of the adjacency matrix of a directed 
weighted graph, that will be refered to as confidence rule (we 
will also refer to the parameters p a ^a' in this way). So the 
confidence rule is a directed weighted graph, whose nodes are 
the opinions in the model (so a model with M opinions would 
have a confidence rule with M nodes) and the arcs represent 
the ways that opinions are allowed to interact. This graph is 
useful as a way of schematizing the opinion interactions and 
as we show in the next sections, it can be used to find the 
properties of the mean-field fixed points. An example of con- 
fidence rule with 4 opinions is given in figure [TJ 




FIG. 1: A confidence rule for 4 opinions. Here pi_>2 = 0.8, 

Pl^3 = 0.5, pi^4 = 0.6, p 2 ^l = 0.25, p 2 ^3 = 1,P3^2 = 

0.17, p4_j.2 = 0.4 and p a ^ a > = otherwise. 



A. The mean-field model 

For the analysis of the mean field case, we consider a vari- 
ant of the Sznajd model, where at each timestep we choose 
q agents at random and if they agree (meaning they are on 
the same state), they attempt to convince r other agents (also 
chosen at random). If the group of q agents has opinion a, 
then each of the targeted r agents is convinced with probabil- 
ity Pa'^a and retains its opinion with probability 1 — p a '^a, 
where a' is the opinion the targeted agent had before the group 
attempted to convince it (and hence it is different for each of 
the r agents). Adding up the probabilities of all possible pro- 
cesses we obtain the mean field equation in the limit of large 
networks: 

ifr = T^llaVa'ivl "Vc'-S-a - Pa^a') , (1) 

a 1 

where r] a is the proportion of sites with opinion a (the de- 
duction of this equation from the underlying Markov chain 
implies that ?/ is actually the expected value of the propor- 
tion) and a time unit corresponds to a Montecarlo timestep 
(MCT), that is, a number of timesteps equal to the number 
of sites in the network. The phase space of this flow is an 
(AI — l)-simplex denoted as SirriM (that is embedded in an 
AI dimensional vector space, in order to make the equations 
more symmetrical), where the vertices correspond to consen- 
sus states and the other states are convex combinations of the 
vertices, with coeficients ?] a : 

where P a is the coordinate of the vertex corresponding to con- 
sensus of opinion a and P is the coordinate in phase space of 
the point representing the state (771, ... , 77^/) (in other words, 
we're using a barycentric coordinate system). 
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The results for the mean-field fixed points can be expressed 
as problems regarding the existence of groups of nodes satis- 
fying certain conditions in the confidence rule and these re- 
sults are the same for all q > 2. We will give here these 
results for a better understanding of the simulations in section 
Hill leaving the mathematical details for later. Because of the 
connection of these results with graph theory, a small glos- 
sary (with examples) will be useful. For the same reason, we 
will interchange freely the notion of a set of opinions with the 
notion of a set of nodes in some graph (like the confidence 
rule). 

B. Graph theory concepts and glossary 

If G is a weighted graph, with adjacency matrix Gi^j 
(that is, the matrix containing the weights of the graph) and 
Gi^j > 0, one can define its directed skeleton Sk^^G) as 
the directed graph with adjacency matrix: 

JO, if Gi_> 3 ■ = 
~\ 1, MGi^j ^0 ■ 

An example of skeleton is given in figure|2] 




FIG. 2: The skeleton for the confidence rule in figureQ] 



• An independent set A is maximal if it contains all the 
nodes in the graph or if the addition of any node not in 
A destroys independence. {3, 4} is a maximal indepen- 
dent set, while {3} is independent but not maximal. 

• If A is a set of nodes from S, then the graph induced by 
A, 6>a is the graph whose set of nodes is A and whose 
connections are the connections between the elements 
of A that existed in S. The graph £{2,3,4} can be found 
in figure [3] 

• The union of 2 graphs G and H, denoted G U H is 
the graph with all the nodes of G and H, but only 
connections that already existed between G and H (in 
short it means referring to 2 unrelated graphs as parts of 
the same graph, without changing anything else). The 
graph 5(1,2,3} U £{4} is shown in figure |4] Also, the 
more familiar concept of component can be defined as 
a graph that is not the union of any smaller parts and is 
also not part of a larger graph with the same property. 

• A graph is strongly connected if we can start at any node 
and get to any other node, respecting the directions of 
the arcs. S, and £{1,2} are strongly connected, but 
<S{i 3} is not because there is no path from 3 to 1 in it. 




FIG. 3: The graph induced in S by the set {2, 3, 4}. 



Let now A be a set of nodes in a directed graph S (typically 
in our problems, A will be a set of opinions and S will be 
the skeleton of the confidence rule). We define the following 
terms (we will use in the examples S equal to the graph in 
figure |2]i : 

• The predecessor set of A, denoted by A_, is the set of 
nodes in S that point to some node in A. For example, 
{3}_ = {1,4}_ = {l,2}and{2}_ = {1,3,4}. 

• Analogously, the successor of A, denoted by A + , is the 
set of nodes in S that are pointed by nodes in A. For 
example, {2,3}+ = {1,2,3}. 

• The complement of A, A is the set of nodes in S that 
are not in A. For example, {2, 3} = {1, 4}. 

• A is an independent set iff S has no connections among 
nodes in A. {1} and {3, 4} are independent sets. Note 
that if A is independent, it follows that A_, A + C A. 



© 



® 



FIG. 4: An example of graph union. The graph Ss\ ^.3} US{4}. 

Finaly, we will denote by Ma the manifold with the states 
where only opinions in A survive: 



Ma = sV G Sim M 



ctGA ) 



(3) 
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C. Mean-field results 

Let TZ be the skeleton of the confidence rule. The results 
for the mean-field fixed points are: 

• Given a fixed point, its stability properties depend only 
on which opinions survive in it and on the skeleton of 
the confidence rule. 

• There exists a fixed point, where all opinions survive, iff 
TZ is a union of strongly connected graphs. Moreover, if 
TZ itself is strongly connected, this point is unique and 
an unstable node (the only exception is the case with 1 
opinion, when the fixed point is the only point in the 
phase space). This is equivalent to the more intuitive 
statement that there exists a fixed point where all opi- 
nions coexist iff we cannot build a set A of opinions that 
can convince opinions in A, but cannot be convinced by 
any opinion in A. 

• The results concerning only opinions in a set A (the 
fixed points and the stabilities inside Ma) can be found 
using the model defined by the confidence rule TZa- (in 
other words, removing opinions from the model leaves 
us with a model with a different confidence rule, that is 
valid inside of Ma)- 

• If 1Z can be split in 2 independent models, that is, 
TZ = TZa U Then for all ff G Ma and V g 
fixed points of the models with rules TZa and TZ-r re- 
spectively, all the points aff+ (1 — a)v with < a < 1 
are fixed points of the model with rule TZ. Moreover, 
the number of unstable directions and stable directions 
is the respective sum of the numbers for ff and V (when 
considering only directions inside Ma and Mj). The 
same thing is true for the neutral directions along which 
there is movement, but fails for the ones with no move- 
ment (we must add ff — v as an extra direction in this 
case). 

• A fixed point where only opinions in A survive is attrac- 
tive iff A_ = A. This also implies that A is a maximal 
independent set (see appendix l A 1 \ and hence that Ma 
is an attractor. 

These results have some interesting consequences and in- 
terpretations, that should be kept in mind when analysing the 
simulation results. 

• The mean-field has no stable situations where 2 inter- 
acting opinions coexist. This means that all possible 
(static) attractors are of the form Ma, where A is a 
maximal independent set. 

• The requirement that A be maximal for Ma to be an 
attractor allows the existence of attractors with surviv- 
ing opinions that do not convince any opinions at all. 

• The condition A_ = A implies that it is possible 
to build confidence rules that have no such attractors. 



These rules display heteroclinic cycles (see appendix 
IA 2b . which cause oscilations with diverging period and 
are heavily affected by finite size effects during simula- 
tions. 

• If every opinion can convince any other (that is 
Pa^a' 7^ for all a ^ a'), then the consensus states 
are the only attractors. 

• In situations where part of the confidence rule can be 
broken in different components, by the removal of some 
opinions, there are manifolds where all points are fixed 
points, and these manifolds can be analised by putting 
together the analysis of each of the components. 



HI. SIMULATION RESULTS AND EXAMPLES 

For our simulations, we used Barabasi-Albert networks 
with 10 5 sites and minimal connectivity equal to 5 (we used 
different networks, but always with these same parameters). 

In order to compare trajectories obtained by simulations 
with trajectories obtained by integrating equations (Q3 we re- 
call that r\ a is the expected value of the proportion of sites 
with opinion a. Because of this and in order to reduce noise, 
we take averages over many simulations (one can also re- 
duce noise by choosing a larger network size). More im- 
portantly, if the initial condition for the mean-field equations 
is (771 , . . . : t)m), then this means that for the corresponding 
simulations, each site must have its opinion chosen at random 
with probability 77^ for opinion a. 

The simulations we will do will be aimed at giving exam- 
ples of the mean-field results from section UTAl some of their 
counter-intuitive aspects and some divergences between the 
simulations and the mean-field. 



A. Attractors and stability 

To illustrate the results about the stability properties of the 
fixed points, consider the rule TZ, depicted in figure [5] (actu- 
aly, a family of confidence rules). The maximal independent 
sets are A = {1, 2},{1, 5}, {3} and {4}, but only {1, 2} and 
{1, 5} obey A_ = A, meaning that the only stationary at- 
tractors are M{ifi} and M{i^}. After a transient we see one 
of 2 situations, the only surviving opinions will be 1 and 2 or 
they will be 1 and 5. We can see this from the time series of 
Vi + V2 + V5 (it tends to 1) and 772 . 775 (it tends to 0, although 
with a longer transient) (figures [6(a) | and |6(b)) . 

The other fixed points can be found by looking at the 
other induced graphs that are unions of strongly connected 
graphs. They are TZ^y, TZ^}, ^-{3,4} and 7\L{i, 2 ,5}- Note 
that 72-/12,5} = ^-{1} U ^{2,5} and that both components are 
strongly connected. This means that we will actualy have a 
line of fixed points connecting some point in the edge P2 P5 
to the vertex Pi. The stability properties of these points are 
in table U A projection of the phase space (where the weights 
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FIG. 5: The skeleton 1Z, of a confidence rule where 2 of the 
4 maximal independent sets generate attractors. These static 
attractors can all be obtained by solving A_ = A in the rule, 
as pointed in section UTAl and are M.{\^} and Ai^i^y. 
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TABLE I: The fixed points of the rule in figure [5]that are not 
in attractors, denoted by the opinions that survive in them (A). 
The line of fixed points connecting the edge P2P5 to the vertex 
Pi is denoted by 1x2,5. For each fixed point, we list the 
number of unstable, stable and neutral directions (u, s and 
n respectively). The relation of these numbers with the sets 
A, A_, A + and the number of components induced by A is 
given in section HV Al 
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FIG. 6: Time series for the rule in figure with weights ei- 
ther or 1, depicting the attractors. Note that in figure |6(a)1 
the ending value is 1 (meaning that opinions 3 and 4 do not 
survive). In figure [6(b)] the ending value is 0, showing that 
opinion 2 and 5 don't survive at the same time. 







FIG. 7: (color online) Phase space projection, depicting the 
structures described in table [Q and the attractors. On the top 
right we see a reordering of the skeleton of the rule making the 
independent sets more evident. The cyan dashed line shows 
the location of the saddle points (1 x 2,5 ; in table HJ), the blue 
shaded trajectories are passing near the point (4), the red ones 
near the point (3) and the green ones near the point (3,4). For 
all these trajectories and the gray ones, lighter shades indicate 
the beginning of the trajectories and darker shades indicate 
their ending. We can see then the trajectories going to the at- 
tractors Aii 1,2} and 7W{i,5}, with some being at first attracted 
by the saddles in (1 x 2,5) before being repelled. We can also 
see the predicted stable direction for the fixed points (3) and 
(4) and the fact that (3,4) is an unstable node. 



in the confidence rule where taken as or 1), showing the at- 
tractors and the features described in this table can be found 
in figure [7] 



B. Surviving inert opinions 

Next, we consider two examples in which we have opinions 
that survive in an attractor, but don't convince any other opin- 
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ion (we'll call them inert). Consider the rules IZi and IZ2 
given in figure[8] In TZ\, ^{1,3} is an attractor, even though 
opinion 1 is inert (can't convince any of the others). In 7Z2, 
■M{4,5}> -M{i,2} and A^{2,3} are attractors, even though opin- 
ion 2 is inert. 




(a) V.! (b) n 2 



FIG. 8: The skeleton of two confidence rules IZi and IZ2, such 
that inert opinions are able to survive in the stationary state. 

We now check that this effect is present in the simulations. 
Time series for the models with confidence rules 1Z\ and IZ2 
(once again, the weights are taken as or 1) are given in fig- 
ures^ to [9(d)] 




TMCs 

(a) {Tii ) 771 + r}3 

FIG. 9: Time series for the rules 1Z\ and 1Z2- 



C. Rules without stationary attractors 

Consider a rule in which all opinions interact (that is, for 
all pair of distinct opinions a and a' either p a ^ a ' ^ or 
Pa' -K? 7^ 0), but such that every opinion a has at least one 
different opinion a' that it can't convince. The independent 
sets of such rule are all unitary, but we have imposed that a' ^ 
{<r}_ and (7 7^ a' , so there are no solutions to A_ = A for 
this rule and hence it has no stationary attractors (it is possible 
to build other types of examples as well). An example for 4 
opinions is given in figure [TOl 

In the appendix I A 21 we prove that if a graph has no so- 
lutions to A_ = A, then it has at least one directed cycle, 
where no edge is doubly connected. In the phase space, these 
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FIG. 9: (cont.) 
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FIG. 9: (cont.) 
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FIG. 9: * 

Time series for the rules TZi and IZ2 in figure[8]with weights either 
or 1. Graphs [9(a)1 and [9(c)1 depict time series containing the full 
attractor (the ending value is either or 1 depending on the attractor 
reached), while graphs [9(b)] and [9(d)1 focus in the surviving inert 
opinion (which always decays, but can reach a non-zero stationary 
value). 



cycles manifest themselves as heteroclinic cycles. These cy- 





FIG. 10: A rule that has no attractors 



FIG. 1 1 : A rule particulary prone to display long transients . 



cles will always be polygonal curves connecting the vertices 
of the simplex that correspond to the nodes the cycle in the 
graph goes through. Moreover, as the cycle goes through the 
nodes in A, it means that A induces in the confidence rule 
a strongly connected graph with one component and as such, 
there exists an unstable fixed point where all the opinions in 
A coexist (in the example of figure [10] the cycle would be 
1— > 2 — > 3 — > A — »• 1, corresponding to the heteroclinic cy- 
cle P1P2P3P4P1 and to an unstable fixed point where all opi- 
nions coexist). Hence, the heteroclinic cycle is fully contained 
in the border of M.a and there is a fixed point in the bulk of 
Ma that leads the trajectories to its border. The typical re- 
sult is that as time goes by, the trajectories get closer to one of 
the cycles, which causes oscillations with a diverging period 
(as they pass each time closer to the consensus states, that are 
fixed points). In simulations, eventually a random fluctuation 
puts the system in a state where one of the opinions in the 
cycle gets extinct, leading the system to a stationary state. 



D. Long transients and stationary states 

In many simulations, there are situations in which the tra- 
jectories get stuck for long times in states that are not attrac- 
tors. In some of these cases, the simulation got to a statio- 
nary state where there are no active connections between the 
agents (that is, a connection between a pair of agreeing sites 
and a neighbour that they can convince, according to the con- 
fidence rule). In other cases there are active connections, but 
some opinions appear in negligible amounts and the set A of 
opinions that are not negligible forms an independent set, but 
not a solution to A_ = A. In the latter cases, the fixed points 
in Ma are saddle points meaning that (this is shown in detail 
in section HVAl A = (A - A_) n A + ^ and usually, one 
(or more) of the negligible opinions will be able to rise again, 
causing long transients. 

Considering the mean-field equations, if a 6 A (meaning 
that a is negligible), then it evolves according to (see section 
IIV Al for further explanations) 



V<r(t) = 



(4) 



1 VtTot^2 a /^:A r 1< T ' P< T '^' J 

as long as the opinions in A are negligible. This implies that 



the time the trajectories spend close to these saddle points can 
be estimated, considering the time it takes for some of the opi- 
nions in A to duplicate its proportion of sites in the network 
(all the other opinions in A will remain negligible for much 
longer times, see appendix iDt. Solving|4]we get 
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eA V 2r 7<ro J2a'eA Va'oPcr' 



(5) 



We now verify this relation for the integration of the mean- 
field equations and compare these results with the simulations. 
We will use the rule in figureQj] with A = {3, 5}. 

For this choice of confidence rule and opinion set, we can 
approximate equation|4]with 



1 

Via 



1 

'/4 



t. 



(6) 



The graphs for the mean field and the simulations can be 
found in figures [T2(a)| and [12(b)| and show that if a simulation 
doesn't go to a stationary state, then it undergoes a transient 
much faster than what is expected from the mean-field equa- 
tions. 



IV. ANALYTICAL RESULTS 



A. Stability and structure of the fixed points 



Let us recall the mean field equation found in section III Al 
(equation [TJ: 



This equation is a flow where the phase space is the simplex 
defined by 



E 



r) 



1 and ricr > OVcr 



and it is easy to show that all trajectories starting in the sim- 
plex never leave it (the sum is because the term be- 
ing summed in equation Q] is skew-symmetrical and the sign 
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(a) Mean Field 
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lit 

(b) Simulations in Barabasi- Albert networks 

FIG. 12: (color online) Return times for A = {3, 5} (more 
precisely the time the trajectory took since it crossed the sur- 
face 773 + 7/5 > 1 — e, until it crossed the surface 774 > 2e, with 
e = 0.025). We can see that simulations behave very differen- 
tly than the mean-field. Particularly, return times are smaller 
than predicted by equation [6] and the relationship between the 
two variables is not linear. The blue points correspond to the 
measured values and the black line corresponds to the predic- 
tion of equation [6] (that holds only for the integration of the 
mean field equations). 
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of the variables doesn't change because dl °g^/dt 
bounded in the bulk of the phase space). 

Let us define if = (771,772,%, ■ • -,Vm), Vc = Fv(v) and 
F = (F 1 ,F 2 ,F 3 ,...,F M ), implying that rf = F(ff). Let 
us also denote the skeleton of the confidence rule by 1Z. The 
fixed points of equation Q] are given by 




Vcr'Pa 



= J2a> Vl'Po 



(7) 



for each opinion a. If if * is a fixed point (meaning, the repre- 
sented point in phase space is a fixed point), then one can look 
at the behaviour of trajectories close to it, using the spectrum 
of the jacobian of F, evaluated at 77*. The jacobian of F, J 
is given by 



+ r6„ itr i 22(qr)l 1 ifr>>Pa>>^o- - r)l"P^a")- (8) 



For each fixed point if*, we can define two sets of opinions, 
A and f2 = A, respectively the set of opinions that survive 
and that are extinct in 77*. It is also useful to look at the com- 
ponents induced by A: Ai , . . . , Afe. That is, 



n A = \Jn Ai . 

i=l 

The jacobian J, can be simplified when it is evaluated at 
the fixed point and il ^ (we will denote all quantities eval- 
uated at the fixed point if * by adding a * superscript): 

J* a , = -rS^, ri*lPa^» , if a £ fi (9) 

rj" 

implying that the jacobian can be written up to permutations 

as 



J* = 



J& 



(10) 



where is the jacobian restricted to the opinions in A and 
i7q is the one restricted to SI (that is a diagonal matrix). It also 
follows from equation [8] that can be written in a block- 
diagonal form where each block J* is the jacobian restricted 
to the component A,: 



J* 



Jt 

J-2 



M 
M- 

Jn. 



(11) 



so in order to find the spectrum of J*, we need to put together 
the spectrum of J7q (which is trivial to find) with the spectrum 
of each of the blocks J* . This means that all that is left is to 
find the spectrum of the jacobian in a coexistence fixed point, 
for a rule with one component. We show in the appendix [C] 
that, if q > 2, all the eigenvalues of J*^ have positive real part 
in this situation and the corresponding eigenvectors are paral- 
lel to the phase space, with the exception of one eigenvector 
that is not parallel and has eigenvalue (with multiplicity 1). 
As this last eigenvector doesn't take points in the phase space 
to other points in the phase space it is irrelevant to the stability 
analysis. This null eigenvalue can be regarded as an artifact 
of having embedded an (M — 1) -dimensional phase space in 
M dimensions. However, according to our reasoning based 
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in equation[TT] it does mean that if A induces fc components, 
there are fc eigenvalues equal to 0. To understand what these 
null eigenvalues mean, consider the case fi = 0. We define 
fji and Fi as the vectors containing only the coordinates of fj 
and F that are in A,-, that is 



rji m 



and F 



We also define the vectors u, as 



vi 



fji 



,V2 



,v k = 



F-l Fn 



fh 



F, 



For every i, Fi is a function of only fji. Moreover, Fi is 
homogeneous with order q + 1, so if 



»7A, = V) r] a and Q = = 0. (12) 

It follows that all the points given by 



k 



a^Cj such that 2_^ ai = ^ an( ^ ctj > V i (13) 



are fixed points of the model. The set of points defined by 
equation [13] is the convex hull of the points defined by the Q 
(we identify these vectors with points in the phase space using 
equation |2]) and as these vectors are linearly independent, it 
means the convex hull must have k — 1 dimensions. 

Euler's theorem gives us another consequence of the homo- 
geneity of Fi, namely 



Jiffi = (?+!)#• 



(14) 



So in the fixed point we have J*fj* = =>• J*v* = 0, yield- 
ing a base of eigenvectors with eigenvalue 0. These eigen- 
vectors can be reorganized so that they are still linearly inde- 
pendent, but only one of them points outwards from the phase 
space, meaning that said eigenvector is irrelevant to the sta- 
bility analysis (the embedding artifact previously discussed). 
Such a base can be taken as Q — fj*, i = 1, . . . , k— 1, together 
with fj* . Removing the artifact, we have a base that generates 
the flat defined by equationQ~3] As all the points in the flat are 
fixed points, this means that there is no movement along these 
directions for trajectories close to the fixed points. 

Finaly, we need to study the spectrum of J7q. As this is a 
diagonal matrix, we can obtain the eigenvalues directly: 



er'GA 



<o, 



(15) 



for each a £ f2 and the eigenvalue \ a tells us if the trajectories 
are attracted or repelled to the manifold r\ a = 0. 



It follows that if TZ is such that SIC A_, then X a < for 
all cr G and we have a complete picture of the behaviour 
of trajectories close to the fixed point. On the other hand, if 
fi ^ A_, some of these eigenvalues are null and we need to 
examine higher orders. In order to do that, if X a = 0, then 
we must make a Taylor expansion of F a around fj* . As F a is 
given by 



o-'ef! 



then it is easy to see that the lower order term, different from 
yields 



er'GA 



(16) 



meaning that the trajectories are repelled from the manifold 
r\ a = 0, unless a ^ A + . In particular, for q = 2 and r = 1, 
the solution of equation[T6lreads 



Va 



as stated in section HH DI (equation[4l). 

Let A4a be the manifold where only opinions in A survive. 
We then know that if all \ a < 0, the trajectories get attracted 
to A4a, but if X a = and a G A + for any a 6 fi, then the 
trajectories get repelled from Ma- Finaly, suppose that either 
A CT < or a £ A + for all opinions in fi. Let u> be the set of 
opinions such that X a = 0, 

A = max Acr < and = > n a . 



We show in the appendixlTJIthat starting in a sufficiently close 
neighbourhood of = the following inequality holds 



|r/„A < 



Ir/A 



(17) 



and so trajectories are neither attracted to nor repelled from 
M A - 

We can now put all these results together. If fj* is a fixed 
point, such that TZa has k components, then the trajectories in 
a neighbourhood of fj* are such that (remembering that fi = 
A): 

• There are |A| + |(A— A_)nA + | — fc unstable directions. 

• There are A n A_ | stable directions. 

• There are fc — 1 directions along which there is no move- 
ment. 

• There are |(A — A_) n (A — A + )| directions that are 
neither attractive nor repulsive, but along which there is 
movement. 
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It follows from the equations Q] that if A is independent, 
then A^a is composed entirely of fixed points. On the other 
hand, if ff * is attractive, then it has no unstable directions and 
the only neutral directions are those along which there is no 
movement, so if k is the number of components of TZa, then 
|A| > k and it follows that 

|A| + |(A- A_)nA + |-fc = 0^ |(A"-A_)nA+| = 

= k- |A| < =>- |(A- A_)nA+| = => |A| = k. (18) 

|A| = k implies that A is independent. The condition that 
there are no neutral directions yields 

(A-A_)nA+ = and (A — A_)n(A — A + ) = 0. (19) 
This implies that 

((A - A_) n A+) U ((A - A_) n (A - A+)) = o 

(A - A_) n (A+ U (A - A+)) = 0. 

As A is independent, A + C A and so A + U (A - A + ) = A, 
implying 

(A-A_)nA = 0«A-A_ = 0«AcA_. 

As A is independent we have A_ C A A = A_. So, 
as A = A_ always solves equations [18] it follows that M.a 
is attractive iff A is independent and A = A_. As we show 
in appendix I A II A = A_ alone implies that A is maximal 
independent, so to find the attractors of the model it suffices to 
find all the sets A, satisfying A_ = A and the corresponding 
manifolds A^a- 

B. Existence of the fixed points 

In the last section, we analysed the stability properties of 
a fixed point, supposing that it existed (given a fixed point, 
what its stability properties are). Now we check when a fixed 
point, where only the opinions in a given set A coexist, exists. 
Our analysis of the case in which the set A of surviving opi- 
nions induces more than one component in the confidence rule 
shows us that we only need to study the case in which all the 
opinions coexist and the confidence rule has one component. 

In this case, any fixed point where all opinions coexist must 
be an unstable node, if it exists, and so if we had embedded 
our phase space in M — 1 instead of M dimensions (substi- 
tuting T]m by 1 — J2a^M for example), the jacobian J of 
the corresponding flux would be a real matrix that is positive 
definite when evaluated in such a fixed point, implying that 
detjj'*) > 0. It follows that the index of the fixed point is 
1 [ 18] and that we can apply the implicit function theorem in 
this case. 



In the appendix|EJ we use these informations together with 
the Poincare-Hopf theorem to show that if our confidence rule 
is such that the directed skeleton is a complete directed graph 
(a graph where there is a doubly connected edge between any 
two nodes) then there exists exactly one fixed point where all 
opinions coexist. On the other hand, a confidence rule 1Z can 
be regarded as a point in the parameter space and for every 
confidence rule, there exists arbitrarily small neighbourhoods 
of it in this space, containing rules whose directed skeletons 
are complete. So for every confidence rule 1Z, there exists 
a path in the parameter space leading to it, but such that all 
other rules in the path have complete skeleton. Finaly, if our 
rule has a complete skeleton, then we can apply the implicit 
function theorem for its coexistence fixed point, meaning that 
the fixed point changes continuously for continuous changes 
in the parameters (changes in the rule). More importantly, the 
jacobian evaluated in the fixed point changes continuously, as 
well as its eigenvalues. 

Suppose then that we have a rule 1Z in which there are no 
coexistence fixed points. We build a path ending in 7Z, where 
all rules have complete skeleton and we look at the coexis- 
tence fixed point along the path. The limit of the fixed point as 
the rule tends to 1Z must also be a fixed point (as F is changing 
continuously too). Examining the eigenvalues of the jacobian 
evaluated at the fixed point along the path, with the exception 
of the limit fixed point, the real parts of the eigenvalues must 
all be positive. For the limit fixed point, there exists a set O of 
non-surviving opinions and Jq must have real non-positive 
eigenvalues, implying that Jq = (by continuity). Recall- 
ing equation [T3J this implies that if A is the set of surviving 
opinions in the limit fixed point, then A_ n A = 0. 

On the other hand, suppose that we have a rule with one 
component and a set A ^ 0, V (V is the set of all nodes) 
such that A_ n A = 0. As we have only one component, we 
must have A + nA ^ 0, and so if we define r/A = So-gA Vcr, 
then 

meaning tja < in the whole region of the phase space where 
all opinions coexist. This implies that there exists a fixed 
point where all opinions coexist iff there is no set of opinions 
A 7^ 0, V, obeying A_ n A = 0. We show in the appendix 
IA 31 that this is equivalent to saying the graph is strongly con- 
nected. So going back to our previous results about the struc- 
ture of the fixed points, we have that there exists a fixed point 
where only the opinions in A survive iff A induces an union 
of strongly connected graphs (as such a fixed point exists iff it 
exists for each of the components separately). 

A similar argument can be used to prove uniqueness. Sup- 
pose that 1Z is a rule with one component that has coexistence 
fixed points. We can build a path ending in 7Z going only 
through rules with complete skeleton. But the coexistence 
fixed points are unique along this path and as we can apply 
the implicit function theorem, any coexistence fixed point of 
1Z must be a limit fixed point of the rules in the path, implying 
that it is also unique. In the case where A induces k strongly 
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connected graphs, this implies that there exists one and only 
one (k— 1) -dimensional flat where only opinions in A survive. 



V. CONCLUSIONS 

On this work, we expanded our previous results about the 
Sznajd model with general confidence rules (interpreted here 
as biases and prejudices), giving analytical results about the 
existence, structure and stability properties of the fixed points 
in the mean-field case, finding a very rich behaviour. We gave 
simulation results in Barabasi-Albert networks that show ex- 
amples of this mean-field behaviour and showed some of the 
discrepances between the model simulated in these networks 
and the integration of the mean-field equations. 

Even though neither the equations for the fixed points can 
be solved analiticaly, nor can the exact eigenvalues of the Ja- 
cobian be all determined, our dynamical systems approach 
was still able to determine the sign of the real parts of these 
eigenvalues and the higher order behaviours, when these were 
needed. Surprisingly, this analysis showed us that the various 
properties of the fixed points depend only on a few qualitative 
properties of the confidence rule (the directed skeleton). This, 
in turn, allowed us to make a connection between the mean- 
field results and graph theory problems and this connection 
can even be used to study more complex behaviours, like the 
heteroclinic cycles in the phase space that always appear in 
the absence of attractors. 

In regard to the simulations, most of the discrepances with 
the mean-field seem to come from the existence of frozen 
states that don't correspond to mean-field attractors, but that 
can be reached by the model in a network. It is not entirely 
clear if these are purely finite size effects, but their origin sug- 
gests that they should be more common as the number of opi- 
nions increases and that the introduction of a random noise, 
in which opinions change randomly with a given probability, 
should destroy this effect. A curious finding in the confidence 
rule studied in section llllDl is that when simulations got close 
to the frozen states, but managed to get away from them, they 
took much less time than would be expected from the mean- 
field results (it must be stressed that we only investigated this 
behaviour for this confidence rule). 

Given the simple conclusions that were reached and the 
generality of our model (we would like to stress that the mean- 
field results are valid not only for the Sznajd model but for 
the g-voter model with q > 2), we believe that similar ap- 
proaches might be fruitful in other models where asymmet- 
rical interactions exist, way beyond opinion propagation and 
sociophysics, like infection spreading and ecology models. It 
would also be interesting to see if similar connections with 
graph theory problems exist in other models and, if they do, 
how rich they are. 
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Appendix A: Graph theory theorems with applications to our 
model 

1. A = A_ implies maximal independence 

Theorem 1. Let G be a graph and let A be a set of nodes 

in it such that A = A This implies that A is a maximal 

independent set. 

Proof. To see this, suppose that A_ = A but A is not inde- 
pendent, then it follows that there exists a, a' £ A such that 
a £ {cr'} + a' £ {o-}_ =^(r'eAnA_/0^ An A ^ 
0, which is a contradiction. 

So if A_ = A then A is independent. If A = then it is 
trivial that A is maximal. If A ^ 0, take a € A. It follows 
that 

(Au{cr})_ = A_U{a}_ = AUW_ => 
(A U {a}) n (A U {<t})_ = (A U {a}) n (A U = 

= ((A u W) n A) u ((A u W) n {*}_) d 

D (AU{tr})nA = {a} ^ => (Ali{a})n(AU{a}) _ ^ 0, 

which implies that A U {a} is not independent and hence, A 
is maximal. □ 

2. Relation between the absence of attractors and heteroclinic 

cycles 

Theorem 2. Let G be a directed graph such that no set of 
nodes obeys A = A_, then there exists a directed cycle in G 
that doesn 't use any of the doubly linked edges. 

Proof. Suppose that there is no such cycle in G and let G 1 be 
the graph G after removing all the doubly linked edges. By 
hypothesis, G 1 is a directed acyclic graph and so a topological 
ordering in G' is possible. This means that we can define a 
strict partial order in V(G): 

i -< j iff there is a path from j to i in G' . 

We can also restrict this order to a subset f2 of V(G), such 
that !,jefl^i -^q j iff i -< j (note that this is not the same 
thing as saying the path exists in G' n ). Consider now the set 
A built from the following algorithm: 

1. Attribute A <- 0, 5 <- and <- V(G). 

2. If fl equals stop, else let i be a minimal element of 

3. Remove i from Q and add it to the set A. 

4. Remove the elements from the predecessor of i with 
respect to G that are in Q, (i_(G) (~l H), from O and add 
them to H. 



12 



5. go to 2. 

By construction the set 5 obeys A = E C A_ (the predeces- 
sor with respect to G). Moreover, the set A is independent. 
To see this, suppose that at some time during the construction 
of A, there are no connections in G between nodes in A and 
nodes in £1 when we reach step 2 (this is trivially true for the 
starting iteration) and let i be the minimal element of -<q cho- 
sen in this step. As i is minimal, there are no nodes in £1 such 
that j <n i and hence there are no paths from i to any other 
element in £1 in the graph G' and hence i+(G) D £1 contains 
only nodes that are connected to i through doubly connected 
edges, implying i+(G) fl O C i-(G) D fl and so in step 4 
we are transfering all the nodes in Q, that had any connec- 
tion with A after step 3, to the set H. So after an iteration of 
the algorithm there are still no connections between nodes in 
A and £1 when reach step 2 again (and so by induction, this 
holds during the whole contraction of A). But as the nodes are 
added to A from £1 one at a time, adding a new node won't add 
connections between nodes in A, implying that A remains in- 
dependent during its whole construction. On the other hand, 
this implies A_ C A. Recalling that by the construction of 3 
we have A C A_, it follows that A_ = A. □ 

The relevance of this theorem to our problem is that a so- 
lution to A = A_ in the skeleton of the rule is equivalent 
to a static attractor in the phase space and saying the cycle 
(j i — >• (72 — > . . . — > <J\ in this skeleton has no doubly con- 
nected edges is equivalent to saying that the polygonal curve 
Pax P<j 2 ■ ■ ■ Pax i s a heteroclinic cycle, meaning that every rule 
that has no static attractors must have at least one heteroclinic 
cycle. 

3. Necessary and sufficient condition for a graph to be strongly 
connected 

Let G be a graph and A ^ a set nodes. 

Def 1. A is a sink (source) iff it obeys A + nA = (A_PA = 
25). In both cases, A is called minimal if there is no non-empty 
proper subset of it with the same property. 

Def 2. The span of a node i, i spa n is the set of all nodes j 
in G, such that j can be reached from i. The span of a set 
of nodes A is defined as the union of the span of each of its 
nodes. 

Corollary 1. Every span is a sink and if A is a sink, then 
A A S p an . 

Corollary 2. As no arc leaves a sink, if X is a sink and Y C 
X then Y span = Y span (G x ), the span ofY in G X - 

Theorem 3. A sink ( source ) is minimal iff it induces a strongly 
connected graph. 

Proof. Let A be a sink that induces a strongly connected 
graph in G. Suppose by absurd that A in not minimal, then 
there exists V C A, such that V is also a sink and V / 0. 
Let lu 6 A — r. As Ga is strongly connected, then for all 



i,j G A we have i £ j spa n{G&) and hence uj e T span (G A )- 
As T C A and A is a sink in G, it follows that r span (GA) = 
r sp an(G) and as T is also a sink in G we have T span (G) = T. 
But this implies ui 6 T, which is a contradiction and so A 
must be minimal. 

On the other hand, if A is a minimal sink in G and we sup- 
pose by absurd that Ga is not strongly connected, there exists 
i,j e A such that i £ j sp an(GA)- A is a sink and {j} C A, 
so this means i <fc j span {G) C and j span {G) C A span (G) = 
A. But then j spa n(G) C A — {i} and so j sp an(G) is a non- 
empty proper subset of A that is a sink, contradicting the as- 
sumption that A was minimal. Hence, Ga must be strongly 
connected. 

The proof for sources is obtained considering the graph 
G', obtained by switching the orientation of all the arcs of G 
(which transforms sinks in sources and vice-versa, but keeps 
the same induced graphs strongly connected) □ 

The relevance of this to our problem is that when the con- 
fidence rule has only one component, the condition that we 
found for the existence of a coexistence fixed point can be 
rephrased as saying that the set of all nodes is a minimal 
source. This theorem shows then that this is equivalent to say- 
ing the confidence rule is strongly connected, which makes 
more easy to see what the result for many components is. 

Appendix B: Topology and matrix theory theorems 

Theorem 4 (Poincare-Hopf). Let M.bea compact, orientable 
and differentiable manifold and let F be a vector field defined 
in M., such that it has only isolated zeros ( every zero has an 
open neighbourhood in which it is unique). If either A4 has no 
border or if F points outwards (acording to the orientation of 
M.) along all points of the border, then the sum of the indices 
of all the zeros of F in the interior of M. equals the Euler 
characteristic of M. 

Theorem 5 (Gershgorin). Let M £ M n (C) be a square ma- 
trix whose general term is rriij. So if X is an eigenvalue of M, 
then there exists an i such that 

|A - m it i\ < 2J 

Theorem 6 (Levy-Desplanques). Let M E M„(C) be an ir- 
reducible square matrix whose general term is rriij. If 

\irn,i\ > |mj,j-| Mi 

and there exists an i such that 

\mi,i\ > ^2 VmA, 

then det(M) ^ 0. 

Theorem 7. Let M G M n (C) be a symmetrical irreducible 
square matrix whose general term is my. If m^i ^ for 
some i, then for all k such that 1 < k < n, there exists an 
irreducible principal submatrix of M with order k. 
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The next theorem is a strengthening of a theorem found in 

Theorem 8. Let M E M n (C) be a square matrix and let 
its hermitian part be H = ( A/ + A/t )/2. If H is positive 
semidefinite, with the multiplicity of equal to fi, then for all 
T 6 M„ (C), such that T is hermitian positive definite, then 
MT (and TM) is positive semidefinite and the sum of the ge- 
ometric multiplicities of its eigenvalues with null real part is 
smaller or equal than fi. 

Proof. By our hypothesis, the eigenvalues of H are non- 
negative real numbers and its eigenvectors can be arranged as 
an orthonormal basis. We can split this basis in 2 parts, {ui}, 
with the eigenvectors with eigenvalue and {v{\, for the oth- 
ers. Define > 0, the eigenvalue such that Hv t = XiVi and 
define U, the linear span of {iii}. Let a; be a column vector 
and x^ its conjugate transpose, so 

x' = aiiii, x" = PjVj, and x = x' + x" => 
* 3 

x^Hx = (z* + x" f )H(x' + x") = 

= {Hx')\x' + x") + x"\Hx') + x'^Hx" = x'^Hx" = 

Hence, x'Hx > x ^ U . On the other hand 
2Re(x t Mx) = xUlx + (xUlx)* = x\M + M^)x = 
2x^Hx. Let S be a nonsingular matrix and w, a normalized 
eigenvector of S^MS, with eigenvalue 7. Taking x = Sw, it 
follows 

Re(7) = Refow^w) = 

= Rc(w^sHlSw) = Re(x f Mx) = x^Hx. 

Define W = S^U = {y e W Sy G U}, so Rc( 7 ) > 
w £ W. As the dimension of U is fx, it follows that 
W also has dimension ^1. The sum a of the geometric multi- 
plicities of the eigenvalues of S'MS with null real part is the 
dimension of the linear span, N, of the corresponding eigen- 
vectors. As all these eigenvectors belong to W and W is a 
linear subspace, then it follows that N is a subspace of W and 
hence a < [i. 

All the properties of the spectrum of a matrix (including 
algebraic and geometric multiplicities) are encoded in its Jor- 
dan canonical form, and this form is invariant by similarity 
transformations, so S'MS, MSS^ and SS'M have the same 
spectral properties. This proves the theorem, as any hermi- 
tian positive definite matrix T can be written as SS\ with a 
non-singular S using a Cholesky decomposition. □ 



Appendix C: Spectrum of the jacobian for a coexistence point in 
a rule with only one component 

Let rf * be a coexistence fixed point (that is, all opinions 
survive) in a model with a rule that has only one component 
and at least 2 opinions. We recall that J*fj* = (due to 
homogeneity) and that if 1 = (1, . . . , 1), then \ J* = (this 
follows from the conservation of the sum of the variables). 
Let D be the diagonal matrix whose diagonal terms are the 
coordinates of ff * (in other words D = diag(?f *)), then the 
symmetric matrix A, defined as 

A = J*D + (J*D) T (CI) 
has off-diagonal terms given by 

A^ = r(l - g)(«,;w CT + vWJPo^>) < (C2) 

and each of the rows (columns) of A sum 0. Moreover, as 
the confidence rule has only one component, A is irreducible, 
implying that at least one off-diagonal term in each row is 
different from and hence all the diagonal terms are positive. 
Finaly, we can use these informations to apply Gershgorin's 
theorem (appendixIBl and find that A is positive semidefinite. 

Denote the principal submatrix of A, obtained by removing 
row and column a by A^ 7 '. If A 0y „i 7^ is an off-diagonal 
term from A, both A^> and A^ a ' are such that one of the 
rows has a positive sum. As every row has an off-diagonal 
term different from 0, then all principal submatrices of A with 
order M — 1 have at least one row that has a positive sum. Ap- 
plying Gershgorin's theorem again, we find that all the A^ 
are positive semidefinite. Moreover, as A is irreducible, then 
there exists a such that A^ is also irreducible and by the 
Levy-Desplanques theorem (appendixIBl. A^ must be posi- 
tive definite. Putting everything together, we find that 

det(A) = and ^dct( J 4 (CT) ) > 0, 

a 

meaning that is an eigenvalue of A with algebraic multiplic- 
ity 1 . As A is also hermitian, then all the other eigenvalues 
must be real and positive. 

As A (a) can be written as A (a) = J* {a) £> (<t) + 
^j*{a) ]j{a)y (because D is diagonal), we can apply our 
theorem from appendix [B] to A^ to find that all the J< a ) 
are positive semidefinite and at least one J*\ a > is posi- 
tive definite. As these matrices are real, this implies that 

det{J*^) > 0. Finaly, we can apply this same the- 
orem to A to find that J* is positive semidefinite and the 
sum of the geometric multiplicities of all eigenvalues with null 
real part is at most one, implying there is at most one eigen- 
value with null real part. As J* .ff* = 0, it follows that is 
the only eigenvalue of J* with null real part. Moreover, as 

det( t 7*( <T - ) ) > it has algebraic multiplicity 1. As we 
also have 1J* = 0, then 1 is a left eigenvector with eigen- 
value 0, implying that if v is a right eigenvector with eigen- 
value different from 0, then l.v = 0. 
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Putting these results together, all the eigenvalues have pos- 
itive real part and the corresponding eigenvectors are parallel 
to the phase space, with the exception of the eigenvector if*, 
that is not parallel and has eigenvalue (with multiplicity 1). 



Appendix D: High order stability analysis for fixed points in 
which opinions get extinct 



-r|n - uj\e qM < — ln(r/ w ) < r|Q - u\e M . 

Integrating in time and taking the limit t — > oo gives the in- 
equalities [T7J 



Suppose that we have a fixed point in which only opinions 
in A survive and let ft = A. For each a £ f2, we define A<j as 



\ a = - r ^ ^ 



/ 4 'la 

cr'eA 



Va'Pa^a' < 0. 



r/qA <n.,< 



r/A 



(D2) 



and so trajectories are neither attracted to nor repelled from 
A4a- This ensures that the whole reasoning is consistent, as it 
is always possible to make rj^o sufficiently small, so that the 
hypothesis of small i] u always holds. 



Suppose that either X a < or a ^ A + for all opinions in 
17 and let uj be the set of opinions such that X a = 0, 



Appendix E: Applying the Poincare-Hopf theorem to the case of 
a confidence rule with complete directed skeleton 



A = max An- < and j] u 



2^ ■no 



It follows from the first order analysis we did in section HV Al 
that if d £ H — oj and the initial value of r\ a , r\ aa is sufficiently 
close to 0, then r] a evolves as 



Consider a rule with a skeleton corresponding to a complete 
directed graph and define 



P : 



min{p 



'}• 



rj a (t) = r\ ao e 

as long as all opinions in £1 remain negligible. It follows from 
the mean field equations that 

^ =7 "E E 11<T'n<T'('n<r 1 P<r'-HT-'>l 9 cr7 1 Pa^<T')- (Dl) 

So if i]u is sufficiently close to 0, it evolves as 

rj a ~ r) ao e Xat V a £ VL - cj => 



E E (r]a'oe^' t 7 1 q a Pa^a-vl' e qK ' t Vai 

(T6u er'Gfi — u> 



yielding the following inequalities 



r E E ^'o f 



E E 

crew cr'£f2— W 



In order to apply the Poincare-Hopf theorem, we build a 
family of manifolds that includes the phase space: 



|?7 £ SirriM r fr > 6 V cr| 



These manifolds all satisfy the hypothesis of the theorem and 
the borders of M e are given by the facets jy^ = e (that is, we 
are using M dimensions to define our manifolds, but we are 
embedding them in M — 1 dimensions). The fixed points we 
obtain for the flow in the mean field equation are not isolated 
when we look at the problem in M dimensions (because of the 
homogeneity of the equations), but our results about the Jaco- 
bian show that embedding the phase space in M — 1 dimen- 
sions instead of M is enough to isolate the zeros (this follows 
from applying the implicit function theorem. Another way of 
isolating the zeros would be to add a term in the equation that 
is inside the phase space, but is different from outside, but 
this has the downside of making the hypothesis to be checked 
more complicated). It also follows from the spectrum of this 
jacobian that the indices of any fixed points in the interior of 
any of the manifolds M e would be 1 . 

The last hypothesis to be checked is then that the vector 
field F points outside along the border. In the manifold M e , 
this is equivalent to the following statement: 



E e 9A -'V<^<r^ J2 eK ' tr >° 



F„ 



Tt E VliPa-HT' < 0. 



(El) 



E E e qA %<V»<rJ2 E 



et-^Vr'Po 



< 



E "^'Po 
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As the left hand side is smaller than e q 1 and the right hand 
side is greater than P/iW 9-1 , then it suffices to take 



e < 



Vc«-i) 
M 



(E2) 



in order to get a manifold M e such that we can apply the the- 
orem. 

The Euler characteristic of all of the M e is 1 , meaning that if 



we can apply the theorem there exists exactly one fixed point 
in its interior. Together with equation lE2] this means that there 
exists exactly one coexistence fixed point and it obeys 



7(9-1) 



M 



Va. 
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